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ABSTRACT 

We report the first spatially-resolved observations of the spectroscopic bina- 
ries A Vir and WR 140, which includes the debut of aperture-synthesis imaging 
with the upgraded three-telescope IOTA interferometer. Using IONIC-3, a new 
integrated optics beam combiner capable of precise closure phase measurement, 
short observations were sufficient to extract the angular separation and orienta- 
tion of each binary system and the component brightness ratio. Most notably, 
the underlying binary in the prototypical colliding-wind source WR 140 (WC7 
+ 04/5) was found to have a separation of ~13 milli-arcseconds with a position 
angle consistent with the images of the 2001 dust shell ejection only if the Wolf- 
Rayet star is fainter than the O star at 1.65/im. We also highlight A Vir whose 
peculiar stellar properties of the Am star components will permit direct testing 
of current theories of tidal evolution when the full orbit is determined. 
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1. Introduction 

Only two optical (visible and infrared) interferometers have published results combin- 
ing three or more separated telescopes, the Cambridge Optical Aperture Synthesis Telescope 
(Baldwin et al. 1996) and the Navy Prototype Optical Interferometer (Benson et al. 1997). 
Combining more than two telescopes is crucial in order to measure the closure phase, a phase 
quantity that is uncorrupted by atmospheric turbulence and necessary for image reconstruc- 
tion in "very long baseline interferometry" (Jennison 1958; Monnier 2000); indeed, closure 
phase lies at the heart of so-called "self-calibration" techniques widely used in the radio (e.g., 
Readhead et al. 1980; Cornwell & Wilkinson 1981). 

Here we present the first scientific results of the upgraded 3-telescope IOTA interfer- 
ometer (Traub et al. 2003), with the new integrated-optics combiner IONIC-3 (Berger et al. 
2003) for precise measurements of visibilities and closure phases. IOTA has much greater 
near-infrared sensitivity (demonstrated H mag 7) than the other imaging arrays, allowing 
new observational programs including imaging of young stellar objects, mass determinations 
of low-mass binaries, and mapping photospheric surface structures of red giants. After de- 
scribing our novel instrumentation, we present the first resolved observations for the A-star 
binary A Vir and colliding wind system WR 140 and discuss the significance of our findings. 
Here we alert the astronomical community to the need for radial velocity follow-up of these 
specific sources as well as to the new research opportunities available with the IOTA3 array. 

2. Observations and Analysis 

All the data presented were obtained using the IOTA (Infrared-Optical Telescope Array) 
interferometer (Traub et al. 2003), which is operated by a consortium of institutions, most 
notably the Smithsonian Astrophysical Observatory and the University of Massachusetts at 
Amherst. The 0.45-m telescopes are movable among 17 stations along 2 orthogonal linear 
arms (telescopes A & C can move along the 35-m northeastern arm, while Telescope B 
moves along the 15-m southeastern arm), allowing an aperture 35mxl5m to be synthesized 
(corresponding to an angular resolution of ~5xl2 milliarcseconds at 1.65/xm). Capping a 
long development (e.g., Schloerb 1990), the third telescope ("C") was fully incorporated on 
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UT 2002 February 23 and now allows visibilities on three baselines and one closure phase 
to be measured simultaneously. The initial observations used a prototype integrated optical 
(10) combiner (e.g., Kern et al. 1997; Berger et al. 2001) which operated only with polarized 
light. The upgraded combiner IONIC-3 corrected this deficiency in 2002 November, and all 
observations presented here were obtained in unpolarized light. 

We only just introduce IONIC-3 (Berger et al. 2003) here; an engineering paper with 
detailed description of the IO component and its performance is in preparation (Berger et 
al. 2004). Light from each telescope is focused into a dedicated single-mode fiber, and the 
three fibers are aligned using a silicon v-groove array mated to the planar waveguides on the 
IO device. The optical circuit acts to split the light from each telescope before recombining 
each telescope pair (AB,BC,AC) at 3 IO couplers. This "pair-wise" combination scheme 
leads to 6-interferometric channels (two for each baseline), where the interference fringes are 
temporally-modulated by scanning piezo mirrors placed in two of the telescope beams; a 
sensitive HgCdTe (Rockwell PICNIC) array (Pedretti et al. 2004a, submitted) detects the 
signals, and data are recorded by a vme-based control system linked to a Sun workstation. 

We followed established observing procedures, as outlined in previous IOTA publica- 
tions (Coude du Foresto et al. 1997; Millan-Gabet et al. 2001): observations of target and 
nearby calibrator stars were interspersed to calibrate slowly varying system visibilities and 
instrumental closure phases. The science targets presented here were observed using a stan- 
dard H band filter (Ao = 1.65/im, AA = 0.30/im); Table 1 contains other observing details. 
Figure 1 shows the (u,v)- cover age obtained during each epoch for each science target. 

Reduction of the V 2 data was carried out using an Interactive Data Language (IDL) 
implementation of the method outlined by Coude du Foresto et al. (1997). In short, we 
measured the power spectrum of each interferogram (proportional to the target V 2 ), after 
correcting for intensity fluctuations and subtracting out bias terms from read noise, residual 
intensity fluctuations, and photon noise (e.g., Perrin 2003). To assure good closure phase 
measurement, we required that interferograms be detected on > 2 baselines, a condition 
nearly always maintained by a realtime fringe packet "tracker" (Pedretti et al. 2004b, in 
preparation). Lastly, the data pipeline applied a correction for the variable flux ratios for each 
baseline by using a transfer matrix (based on single-beam flux measurements; e.g., Coude 
du Foresto et al. 1997; Monnier 2001). We have studied the absolute calibration accuracy 
by observing single stars of known size and found a systematic error of AV 2 = 0.10 (e.g., 5% 
error in visibility for unresolved sources) under the worst observing conditions; this error 
has been combined (in quadrature) with the statistical error for fitting purposes. Improved 
analysis methods should achieve <1% errors eventually, since we are using single- mode fibers. 

We followed the method of Baldwin et al. (1996) for calculating the complex triple 
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amplitude in deriving the closure phase, explicitly guaranteeing "fringe frequency closure" 
{vab + vbc + v ca = 0). While photon-noise bias correction is not necessary for the IONIC-3 
( "pair- wise" ) combiner, the instrumental closure phase must be determined empirically. The 
miniature dimensions of the 10 component have minimized drifts to less than 1 degree over 
many hours; chromaticity effects, however, limit our absolute precision when the calibrator 
and source are not of the same spectral type due to different "effective wavelengths." En- 
gineering tests indicate that the instrumental closure phase (<&cp) varies systematically by 
1.4 ± 0.3 degrees between a hot star (B8) and a cool star (M3) when using the broadband 
H filter. While we minimize these errors by using calibrators matched to the spectral type 
of the target, a conservative systematic error of A$ CP = 0.75 degs was adopted here to ap- 
proximate our expected residual level of miscalibration. Lastly, the sign of our closure phase 
(J&ab + &bc + &ca for baseline triangle connecting telescopes A— »B— >C— >A) was calibrated 
using the well-known binary stars Capella and Matar (Hummel et al. 1994, 1998). 



3. Results and Discussion 

With calibrated V 2 and closure phases for each target, a static binary star model (sepa- 
ration, position angle, brightness ratio) was fitted to the data for each given epoch, assuming 
the calibration errors given in the previous section and ignoring bandwidth-smearing effects. 
Table 2 contains the results of our fits based on a maximum-likelihood approach giving equal 

weights to the visibility data and closure phases. Explicitly, the probability of a given set 

x 2 2 

of parameters was approximated by: Probability oc exp( — — ^ — r ■ ■ Xcv , i \ )- This 

1 11 J mimmum(x' 2 ) mmimum(Xpp)' 

normalization is equivalent to increasing the data errors to force the x 2 /DOF= 1 ? and con- 
servatively estimates errors in the context of an incomplete error model. The most probable 
values with 1-sigma errors bars (68.4% confidence interval) were calculated for each param- 
eter (separation, position angle, brightness ratio) by integrating over the other variables, a 
standard practice referred to as marginalizing over these variables. 

Based on their brightness and color temperatures, the component stars of both bina- 
ries are expected to be unresolved here (^0.5 mas); this assumption, as well as neglecting 
bandwidth-smearing, may not remain valid if the data precision is improved significantly. 
The V 2 and closure phases for a subset of the data are shown in Figure 2 along with the 
best-fitting model prediction. We note that A Vir, which has a ~ 207 day period, is not 
completely static during each observing "epoch." With complete orbital phase coverage, this 
deficiency will be corrected by solving for the full orbital elements by fitting directly to the 
visibilities, closure phases, and radial velocities, a procedure outlined in detail by Hummel 
et al. (1998) and Boden et al. (1999). 
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The radio community developed interferometric imaging techniques (e.g., hybrid map- 
ping) capable of reobtaining lost individual phases (e.g., Haniff et al. 1987). As demonstrated 
by Baldwin et al. (1996), these methods can be applied successfully to optical interferometers 
as well. We present the first aperture synthesis images from IOTA3 in Figure 3, depicting 
A Vir for each epoch. The images are constructed iteratively starting with zero initial phases 
and using the measured closure phases for self-calibration. CLEAN deconvolution was ap- 
plied using an elliptical Gaussian CLEAN beam. Full details of the image-making process 
for this source (as well as the binary Capella) can be found in Kraus (2003). The presented 
maps are consistent with the binary parameters found through direct fitting, and demon- 
strate the imaging capability of the IOTA3+IONIC-3 system, currently the most sensitive 
infrared imaging interferometer in the world. 

With this new imaging capability, we begin to study binary stars which have never been 
resolved before and which represent unique stellar systems whose properties can shed new 
light on rare and/or short-lived phases of stellar evolution. 

A Vir is a well-known spectroscopic binary (period 207 days, Stickland 1976) with similar 
components of spectral type early Am (both are metallic-lined, Abt 1961). The stars exhibit 
the curious property that one is sharp-lined and the other broad-lined, which makes them 
easily distinguishable spectroscopically. IOTA3 observations will allow a determination of 
the absolute properties of the components once the full orbit is established. Fundamental 
stellar parameters, such as the system distance (independent of Hipparcos) and the individual 
masses and absolute near-IR luminosities, can be determined through a combination of the 
radial velocity orbit, astro metric orbit, and brightness ratio. Further, these properties can 
be compared directly against models of stellar evolution to derive the age of the system. The 
very different projected rotational velocities of the components is of special interest for these 
Am stars, enabling unique tests of current theories of tidal evolution (Tassoul & Tassoul 
1997; Zahn & Bouchet 1989, and references therein). 

WR 140 is a prototypical colliding-wind binary source consisting of a Wolf-Rayet (WC7) 
and an 04-5 star in a highly-eccentric ~7.9 yr orbit, and has been extensively studied in the 
radio, infrared, optical, ultra-violet, and X-rays (e.g., Moffat et al. 1987; Williams et al. 1990; 
White & Becker 1995). Notably, infrared emission appears and fades in conjunction with 
periastron passage, explained by transient dust formation in the dense shock- compressed gas 
at the wind-wind interface (Williams et al. 1990; Usov 1991; Tuthill et al. 1999). This process 
was captured in a series of high- resolution near-infrared images by Monnier et al. (2002), 
revealing an eastward moving arc of dust consistent with the O-star being roughly east of 
the WR star at periastron. The O-star is expected to be on (nearly) the opposite side of the 
WR star now as compared to periastron, which can be reconciled with the position angle 
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found here (PA 152°) only if the WR star is fainter than the O-star at 1.65//m (as expected 
by some workers; e.g. Williams et al. 1990). We intend to combine future high-precision 
interferometer observations with new radial velocity data (e.g., Marchenko et al. 2003) to 
derive component masses and a distance estimate, fundamental parameters for these systems 
needed for definitive interpretation of data across the electromagnetic spectrum. 

4. Conclusions 

We have reported first scientific results from the IOTA3 interferometer: newly-resolved 
separation, position angle, and flux ratio information for double-lined spectroscopic binary 
stars A Vir and WR 140. Using a preliminary data reduction pipeline, calibration is sufficient 
for model-fitting and aperture-synthesis imaging using hybrid mapping techniques. New 
observations are being pursued in order to derive precise component masses and distance 
estimates by combining interferometric and radial velocity data. All the visibility and closure 
phase data here have been converted to the new FITS format for Optical Interferometry data 
(OI-FITS) 1 and are available upon request. 
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Table 1. Observing Log for A Vir and WR 140 



Date 


Interferometer 


Calibrator Names 


(UT) 


Configuration ^ 


(Adopted UD Diameter) 




A Vir (A1V + A); H mag = 4.28±0.21 


2003 Feb 16,17 


A35-B05-C10 


HD 126035 (G7III, 0.78±0.24 mas) 






HIP 71957 (F2III, 1.20±0.22 mas) 


2003 Feb 20,22,23 


A25-B15-C10 


HD 126035 


2003 Mar 21 


A35-B07-C25 


HD 126035 


2003 Mar 22 


A35-B07-C10 


HD 126035; HD 158352 (A8V, 0.44±0.10) 


2003 Mar 23-24 


A35-B15-C10 


HD 126035; HD 158352 


2003 Jun 12-17 


A35-B15-C10 


HD 126035 



WR 140 (WC7 + 04/5); H mag = 5.429±0.023 
2003 Jun 17 A35-B15-C10 HD 192985 (F5IV, 0.46±0.10 mas) 

HD 193631 (K0, 0.44±0.10 mas) 



Notes: H-band photometry from 2MASS All-Sky Release; calibrator diameters estimated 
using getCal software (http://msc.caltech.edu). 

Configuration refers to the location of telescopes A,B,C on the NE, SE and NE arm 
respectively; see §2 for more details. 



Table 2. Binary Parameters for A Vir and WR 140 at 1.65/im 



Source 


Epoch 


Separation 


Position Angled Flux 


xVdofw 




(U.T.) 


(milliarcseconds) 


(degrees) Ratio 


V 2 CP 



A Vir 2003 Feb 16-23 
2003 Mar 21-24 
2003 Jun 12-17 

WR 140 2003 Jun 17 



18.2 ±0.6 
17.8 ±0.5 
19.5 ±0.6 
12.9lg;t 



200 ±3 
184 ±3 
22 ±3 



151.7; 



-1.8 
-1.3 



1.69lg;^ 0.45 

1.72 ±0.07 0.33 

2.051°!? °- 31 

i-io±8:8i o-is 



12.8 
2.2 
4.8 
4.7 



a Position angles (East of North) are measured from the bright to the faint component. 

b The x 2 P er DOF (degree of freedom), the reduced x 2 , is reported for the model fit to the 
V 2 and closure phase (CP) measurements separately. 



30 20 10 -10 -20 30 20 10 -10 -20 
U Coordinate (m) U Coordinate (m) 



Fig. 1. — (u,v)-plane coverage of science binaries. Each diamond symbol represents a single 
observation with the IOTA3 interferometer, corresponding to approximately 100-200 inter- 
ferogram scans. The overlaid image in each panel represents the V 2 of the best fit binary 
model for each epoch discussed in text. 
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Fig. 2. — (top panel): This figure shows example binary model fits for A Vir on 
UT 2003 March 23. The top three panels show the V 2 for baselines connecting telescopes 
A-B, B-C, and A-C, respectively. The bottom panel shows the calibrated closure phase for 
the A-B-C triangle, (bottom panel): Same as top panel for UT 2003 June 17 observations 
of WR 140. These curves represent the best-fit model parameters as described in §3. For 
all panels, the error bars shown represent only the statistical error and do not include the 
systematic error discussed in §2. 
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Fig. 3. — First aperture synthesis images from the upgraded IOTA 3-telescope interferom- 
eter, depicting A Vir based on 1.65//m observations at three separate epochs. The binary 
components are well-resolved from each other but individually are unresolved; orbital motion 
is clearly evident, showing nearly a ~180° rotation between the first and last epoch. Each 
panel includes a contour map showing 5% intervals in surface brightness scaled to the peak, 
as well as the the CLEAN restoring beam as inset. North is up and East is left. 



